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Based  on  the  fluxes,  energy  spectra,  and  pitch-angle  distributions  of  the  part- 
icles, ion-pair  production  profiles  were  calculated  and  electron  density  pro- 
files subsequently  obtained  with  application  of  known  effective  electron  loss 
rates.  Prom  the  coordinated  data,  correlative  studies  between  signal  strengths 
and  electron  density  profiles  were  performed  with  use  of  an  ELF  wave -guide -mode 
computer  program  developed  at  the  Naval  Ocean  Systems  Center.  With  the  inclu- 
sion of  data  from  several  stations,  the  transmission  signal  strength  could  be 
studied  simultaneously  under  a variety  of  ionospheric  conditions.  The  wave- 
guide computer  program  calculations  indicated  that  energetic  particle  precip- 
itation events  may  cause  a signal  enhancement  or  degradation,  depending  upon 
where  the  disturbance  is  occurring  along  the  path.  The  consistent  occurrence 
during  electron  precipitation  events  of  signal  strength  anomalies  in  either 
direction  is  supported  by  the  coordinated  satellite/ELF  station  measurements. 

A 

The  following  major  conclusions  have  followed  from  this  study  effort: 

• From  both  the  satellite  and  ELF  station  measurements  it  has  been  found 
that  direct  particle  precipitation  into  the  atmosphere  can  cause  ELF 
transmission  anomalies.  The  signal  strengths  may  be  either  attenuated 
or  enhanced  depending  upon  the  geometry  and  details  of  the  ion  pair 
production  profiles. 

• Sensitivity  studies  have  been  made  to  assess  the  dependence  of  the  ELF 
signal  strengths  on  such  parameters  as  the  electron  and  ion  density 
profiles  and  their  distribution  along  the  propagation  path. 

• The  signal  strengths  tend  to  decrease  with  increasing  electron 
density  at  altitudes  of  ~ 60  km  or  lower. 

• The  effect  of  a given  ionization  profile  depends  strongly  on  its 
location  along  the  propagation  path. 

• The  ELF  signal  strengths  are  most  sensitive  to  positive  ion  density 
profiles  at  altitudes  of  ~ 45  km  or  lower. 

• The  ELF  signal  strengths  are  very  sensitive  to  sporadic  E-layers, 
with  the  altitude  of  the  ledge  being  a very  critical  parameter. 

• The  measured  and  predicted  effects  of  energetic  electron  precipitation 
events  can  provide  a readily  available  verification  of  the  effects  on 
ELF  transmission  of  more  rarely  occurring  and  possibly  more  severe 
phenomena  such  as  polar  cap  absorption  (PCA)  events. 

• Variations  in  the  nighttime  ELF  signal  strengths  on  a fine  time  scale 
are  observed  which  may  be  due  entirely  to  electron  precipitation  but 
cannot  be  accounted  for  quantitatively  due  to  present  limitations  in 
the  measurements  and  computational  techniques. 


• The  geometry  for  the  effect  of  electron  precipitation  on  nighttime  ELF 
transmission  is  very  complex  and  as  a result  the  following  recommendations 
are  made  for  future  investigations: 
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• New  techniques  for  mapping  electron  precipitation  profiles 
simultaneously  over  a broad  spatial  region  should  be  used 
in  a coordinated  measurement  program. 

• Existing  ELF  waveguide-mode  computer  programs  should  be  modified 
to  include  treatment  of  variations  in  the  electron  and  ion  density 
profiles  along  a direction  perpendicular  to  the  propagation  path. 
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PREFACE 


The  transmissions  of  extremely-low-frequency  (ELF)  waves  from  the  U.  S. 

Navy  operated  Wisconsin  Test  Facility  (UTF)  to  a receiving  site  in  Connecti- 
cut are  known  to  experience  significant  signal  strength  reductions  of  up  to 
3 dB  on  approximately  60  to  80  nights  per  year. . Transmissions  to  other 
stations  in  Maryland,  Greenland,  Norway  and  Italy  display  similar  anomalies. 

These  signal  strength  reductions,  which  may  result  in  degraded  system  performance, 
have  been  speculated  to  be  caused  by  the  precipitation  of  electrons  from  the 
radiation  belts  resulting  in  enhanced  ionization  in  the  D-  and  E-region  ion- 
osphere. In  the  past,  attempts  to  correlate  the  anomalous  signal  strength 
behavior  with  geomagnetic  indices  indicative  of  disturbed  magnetospheric  con- 
ditions have  met  with  limited  success.  These  inconclusive  studies  clearly 
pointed  out  the  need  to  achieve  a convincing  cause-and-effect  relationship 
by  comparing  the  ELF  signal  strengths  directly  with  particle  precipitation. 

The  opportunity  for  performing  such  a correlative  study  was  provided  by  data 
acquired  with  scientific  payloads  on  the  low-altitude  polar  orbiting  satellites 
1971-089A  and  1972-076B  developed  by  the  Space  Sciences  Laboratory  of  the 
Lockheed  Palo  Alto  Research  Laboratory  (LPARL)  for  the  Office  of  Naval  Research, 
the  Defense  Nuclear  Agency  and  the  Defense  Advanced  Research  Projects  Agency. 

The  latter  satellite  is  still  operational  and  capable  of  acquiring  coordinated 
data.  In  the  course  of  achieving  the  prime  goals  of  the  two  different 
satellite  experiments,  a large  body  of  scientific  data  had  been  obtained  on 
the  precipitation  of  energetic  particles  into  the  atmosphere.  To  investigate 
the  relationship  between  the  ELF  signal  strengths  and  the  precipitating 
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electron  fluxes  using  these  satellite  data,  the  Office  of  Naval  Research 
supported  an  effort  by  the  LPARL  group  to  analyze  the  satellite  data  in  the 
vicinity  of  and  at  the  times  of  measured  signal  strength  reductions  in  the 
WTF-to-Connecticut  link.  Under  the  initial  phase  of  this  program  a successful 
correlation  was  made  between  direct  particle  precipitation  and  anomalous  ELF 
signal  strength  (imhof  et  al.,  1976). 

On  all  of  the  11  coordinated  cases  studied  in  the  first  phase  of  this 
investigation  significant  fluxes  of  energetic  electrons  having  substantially 
different  intensities  and  spectral  shapes  were  observed  to  be  precipitating 
into  the  atmos-  Resulting  ion -product ion  rates  in  the  upper  D-region 

(70-100  km)  xated  to  be  often  comparable  to  those  produced  by  in- 

tense solar-particle-events  (SJE) . Corresponding  enhanced  electron  density 
profiles  were  at  times  in  excess  of  lCre  /cmJ  at  altitudes  above  ~ 80  km.  A 
positive  qualitative  correlation  was  therefore  established  for  the  first  time 
between  anomalously  low  nighttime  ELF  signal  levels  and  the  precipitation  of 
significant  fluxes  of  electrons  from  the  radiation  belts,  but  a detailed 
quantitative  correlation  between  the  severity  of  the  anomalies  and  particle 
inputs  remained  to  be  performed. 

From  the  previous  effort  it  became  clear  that  for  a given  transmission 
site  the  signal  amplitudes  received  at  various  stations  around  the  world  may 
vary  with  time  in  an  appreciably  different  manner.  If  the  fluctuations  are 
associated  with  energetic  particle  precipitation  then  the  station-to-station 
variations  may  reflect  latitude  and  longitude  changes  in  the  total  precipi- 
tating particles  and  energy  spectra.  Stimulated  by  these  considerations  and 
the  previous  findings  a special  coordinated  exercise  was  conducted  in  March- 
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April  1976  between  LPARL  satellite  measurements  of  the  precipitating  particles 
and  ELF  transmission  between  the  U.  S.  Navy  Wisconsin  Test  Facility  and 
receiving  stations  in  Maryland,  Greenland,  Norway  and  Italy  coordinated  by 
Dr.  John  Davis  of  the  Naval  Research  Laboratory.  Also  a large  body  of  satellite 
data  has  been  acquired  during  1974  and  1975  in  the  general  vicinity  of  the  ELF 
paths  to  receivers  at  the  various  stations  during  times  of  operation,  and 
analysis  of  these  data  would  clearly  be  of  great  value. 

The  purpose  of  this  report  is  to  present  the  findings  of  a study  using 

available  satellite  data  taken  during  nights  when  the  ELF  receivers  were  in 

operation  and  experiencing  both  normal  and  anomalous  signal  strength  behavior. 

For  selected  cases  comparisons  were  made  between  the  measured  signal  strengths 

and  theoretical  results  using  the  ELF  waveguide -mode  computer  program  developed 

at  the  Naval  Ocean  Systems  Center  (NOSC) . This  study  has  shown  for  example, 

that  an  intense  Solar  Particle  Event  (SPE)  such  as  the  one  measured  by  LPARL 

at  1144  UT  on  4 August  1972  would  decrease  the  ELF  signal  strengths  from  WTF 

to  Norway  by  2-5  dB  for  a nighttime  path.  Unfortunately,  no  cooperative  ELF 

measurements  were  made  during  this  highly  disturbed  period  that  effected  the 

, o 

entire  polar  cap  down  to  magnetic  latitudes  of  ~ 60  . On  the  other  hand,  for 
various  hypothesized  electron  and  positive  ion  density  profiles  based  on 
variations  of  the  measurements  made  at  1254  UT  5 August  1972,  when  the  solar 
particle  intensities  were  much  lower  than  at  the  peak  on  4 August  1972,  the 
calculated  signal  strengths  are  found  to  decrease  with  increasing  electron 
density  at  60  km,  but  the  signal  levels  may  be  above  the  ambient  levels. 

For  relativistic  electron  precipitation  events  the  coordinated  measurements 
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and  calculations  indicate  that  signal  enhancements  or  decreases  may  result, 
depending  upon  the  details  of  the  density  profiles  and  where  the  disturbance 
is  occurring  along  the  path.  The  experimental  results  and  the  calculations 
of  these  effects  are  discussed  in  this  report. 

Many  individuals  contributed  to  this  successful  program.  Special  ac- 
knowledgments are  extended  to  Mr.  R.  G.  Joiner  of  the  Office  of  Naval  Research 
and  to  Dr.  T.  Quinn,  while  at  the  Office  of  Naval  Research,  for  their  im- 
portant cooperation,  support  and.  direction  under  Contract  N00014-75-C-0954. 

Dr.  J.  R.  Davis  of  the  Naval  Research  Laboratory  generously  provided  all  of 
the  experimental  ELF  data  used  in  this  study.  Mr.  W.  Moler  of  the  Naval 
Ocean  System  Center  (NOSC)  kindly  cooperated  in  providing  the  NOSC  wave 
propagation  code.  We  acknowledge  the  cooperation  of  the  Norwegian  Defence 
Research  Establishment  in  granting  Dr.  T.  Larsen  a leave  of  absence  to  engage 
in  this  ELF  study  activity  while  in  residence  at  LPARL. 
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Section  1 
INTRODUCTION 

In  recent  years  it  has  been  recognized  that  at  least  some  of  the  ob- 
served ELF  propagation  anomalies  in  the  U.  S.  Navy  experimental  transmission 
link  are  due  to  anomalous  ionospheric  ionization  caused  by  the  precipitation 
of  energetic  electrons  and  protons  into  the  earth's  atmosphere  (Davis,  1974; 
Davis  and  Meyers,  1975;  Larsen,  1974).  For  the  first  time  in  1976  a success- 
ful correlation  was  made  between  direct  particle  precipitation  and  anomalous 
nighttime  ELF  signal  strengths  (imhof  et  al.,  1976).  In  the  latter  study, 
energetic  particle  data  were  obtained  with  scientific  payloads  on  the  low- 
altitude  polar-orbiting  satellites  1971-O89A  and  1972-076B.  During  all  11 
coordinated  cases  studied,  significant  fluxes  of  energetic  electrons  having 
substantially  different  intensities  and  spectral  shapes  were  observed  to  be 
precipitating  from  the  radiation  belts  into  the  atmosphere.  Resulting  ion- 
production  rates  in  the  upper  D-region  were  calculated  to  be  often  comparable 
to  those  produced  by  intense  solar-particle-events.  Corresponding  enhanced 

c O 

electron  density  profiles  were  at  times  in  excess  of  lCr  electrons/cm  . A 
positive  qualitative  correlation  was  therefore  stablished,  but  a detailed 
quantitative  correlation  between  the  severity  of  the  anomalies  and  particle 
characteristics  such  as  intensity  and  spectra  was  not  firmly  established 
because  of  the  limited  data  base  used. 

In  this  effort  performed  under  Contract  N00014-75-C-0954  with  the  Office 
of  Naval  Research  the  satellite  data  taken  during  special  coordinated  exer- 
cises in  1975  and  1976  were  analyzed  and  compared  with  measured  ELF  trans- 
mission between  the  U.  S.  Navy  Wisconsin  Test  Facility  and  receiving  stations 
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in  Maryland,  Greenland,  Norway  and  Italy.  The  times  of  the  satellite  passes 
used  in  the  study  are  listed  in  Section  2 along  with  a brief  description  of 
the  payload.  In  Section  3 the  waveguide  mode  program  is  described  along 
with  a discussion  of  the  path  segmentation  used.  Subsequently  the  satellite 
data  and  the  results  of  the  signal  strength  calculations  are  presented  for 
the  March-April  1976  coordinations  (Section  4),  for  other  relativistic  elec- 
tron precipitation  cases  in  1975  (Section  5)5  and  for  Solar  Particle  Events 
(Section  6) . In  Section  7 waveguide  mode  calculations  are  presented  for 
selected  sporadic  E-layer  profiles.  In  Section  8 a short  summary  is  given 
along  with  the  conclusions  of  the  present  study. 
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Section  2 
THE  SATELLITE  DATA 

2.1  The  LPARL  Particle  Payload 

Under  this  program,  direct  correlations  have  been  made  between  ELF 
signal  strength  observations  and  coordinated  measurements  from  a satellite 
of  particles  precipitating  into  the  atmosphere.  The  energetic  particle  data 
were  obtained  by  the  Space  Sciences  Laboratory  at  LPARL  in  two  different  ONR 
sponsored  satellite  programs  having  the  designations  1971-089A  and  1972-076B. 
Each  of  these  satellites  has  the  advantage  of  being  at  low  altitudes 
(750-800  km)  and  polar  orbiting.  Most  of  the  coordinated  data  used  in  the 
analyses  reported  here  were  taken  with  the  satellite  1972-076B  which  is  spin- 
ning and  the  instrumentation  can  therefore  provide  good  pitch  angle  distri- 
bution measurements.  This  satellite  is  still  operational  at  the  time  of  this 
report . 

An  example  of  a coordinated  measurement  is  shown  schematically  in 
Figure  2-1.  The  propagation  path  from  the  transmitter  at  the  Wisconsin  Test 
Facility  (WTF)  to  a receiver  at  Tromso,  Norway  is  illustrated.  At  the  time 
of  interest  the  satellite  passed  through  the  outer  radiation  belt/auroral 
zone  region  where  intense  fluxes  of  energetic  electrons  were  observed  to  be 
precipitating  into  the  atmosphere.  The  fluxes  of  electrons  above  160  keV 
are  shown  schematically.  It  is  important  to  measure  the  particle  fluxes 
simultaneously  with  conduction  of  the  ELF  transmissions  since  the  fluxes  and 
energy  spectra  of  precipitating  electrons  are  known  often  to  be  very  dynamic 
in  time  and  space.  Since  variations  with  longitude  also  occur,  data  from 
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multiple  satellite  crossings  were  used  whenever  available.  In  both  satellite 
experiments  high  data  coverage  was  achieved  on  a worldwide  basis  through  use 
of  tape  recorders.  On  each  of  these  satellites  the  LPARL  payload  contains 
an  array  of  instruments  measuring  electrons,  protons,  and  alpha  particles 
over  a broad  range  of  energies.  A gamma  ray  detector  in  the  1972-076B  pay- 
load  has  also  provided  data  on  widespread  electron  precipitation  events  through 
observation  of  the  bremsstrahlung  x-rays  (imhof  et  al.,  1974a).  These  data 
provide  measurements  of  the  spatially  integrated  electron  fluxes  incident  over 
a large  area,  but  due  to  depletion  in  May  1973  of  the  cryogen  associated  with 
that  instrument  the  data  are  not  applicable  to  most  of  the  event  of  interest 
here . 

Of  particular  concern  here  are  the  data  taken  with  the  energetic  elec- 
tron spectrometers  measuring  the  intensities  and  energy  spectra  of  electrons 
from  ~ 130  keV  to  2.5  MeV  with  256-channel  analysis.  Two  such  identical  elec- 
tron spectrometers  are  located  on  the  oriented  satellite  1971-089A  at  look 
angles  of  90°  and  20°  with  respect  to  the  local  zenith.  On  the  spinning 
satellite,  1972-076B,  a single  spectrometer  provides  detailed  pitch  angle 
distributions  on  the  same  time  scale  as  the  satellite  spin  period,  5 seconds. 

The  Lockheed  group  has  obtained  data  on  precipitating  electrons  from  one  of  these 
two  low  altitude  satellites  continuously  since  1971  (imhof  et  al.,  1973; 

1974b,  1974c). 

In  addition  to  the  energetic  electron  measurements,  on  each  of  the  f-wo 
satellites  the  LPARL  payload  also  provided  measurements  of  the  intensities 
and  energy  spectra  of  lower  energy  electrons  and  of  protons  over  a broad 
range  of  energies.  The  low-energy  (auroral)  particle  complement  on  the 
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1971-089A  spacecraft  consisted  of  34  individual  detectors  measuring  electrons 
in  the  energy  range  from  6 «V  to  40  keV  and  protons  in  the  energy  range  from 
0.6  keV  to  greater  than  *40  keV.  These  detectors  were  arranged  into  six 
instrument  packages  including  three  multiple  particle  analyzers  and  three 
scintillator-photomultiplier  instruments.  A mass  spectrometer  capable  of 
distinguishing  protons,  energetic  helium  ions  and  heavier  ions  over  the 
energy-per-unit-charge  range  0.7  keV  to  12.1  keV  was  also  included.  On  the 
1971-089A  spacecraft  detailed  fluxes  and  spectra  of  energetic  protons 
(1.2  to  100  keV) , alpha  particles  (6.0  to  60  keV)  and  to  a lesser  extent 
electrons,  (O.76  to  > 2 keV) , were  obtained  with  two  spectrometers.  These 
were  essentially  identical  instruments,  one  oriented  to  view  along  the  zenith 
direction  continually  while  the  other  was  oriented  at  90°  to  this  direction. 
The  particle  instruments  in  the  0NR-001  payload  on  the  1971-089A  satellite 
are  summarized  in  Table  2-1. 

The  LPARL  payload  on  the  spinning  satellite  contains  an  array  of 
particle  detectors  in  addition  to  the  Energetic  Electron  Monitor  (EEM)  using 
a plastic-scintillator-photomultiplier  detector  and  a 256-channel  analyzer 
for  fine  energy  resolution.  The  lower  energy  environment  is  measured  with  a 
multi-particle  analyzer  containing  a total  of  nine  channel-multiplier  detec- 
tors. In  this  instrument  three  of  the  detectors  are  devoted  to  protons  from 
1 to  200  keV  and  six  detectors  to  electrons  from  0.06  to  50  keV.  Two  spec- 
trometers measure  the  intensities  and  energy  spectra  of  protons  from  1 to 
100  MeV,  electrons  from  0.8  to  2.1  MeV  and  alpha  particles  from  6 to  *400 
MeV.  A low  energy  proton  detector  (LEp)  measures  protons  over  the  energy 
range  from  0.2  to  1 MeV  in  four  energy  channels.  In  addition,  an 
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omnidirectional  detector  measures  the  integral  electron  flux  above  1 MeV. 
The  instruments  on  the  1972-076B  satellite  are  summarized  in  Table  2-2. 


i 


2.2  Electron  Precipitation  at  Middle  Latitudes 

Energetic  electron  precipitation  is  known  to  occur  on  a worldwide 
basis.  The  precipitation  patterns  are  very  complex,  being  strongly  dependent 
upon  longitude,  latitude  and  geomagnetic  activity.  Many  studies  have  been 
made  of  energetic  electron  precipitation,  but  in  general  concerted  efforts 
have  not  been  made  to  coordinate  the  energetic  electron  observations  with 
the  ionospheric  or  atmospheric  phenomena  of  interest,  such  as  ELF  wave 
propagation  observations,  since  the  fluxes  can  change  on  a very  short  time 
or  small  spatial  scale  it  is  important  to  measure  the  electrons  very  close 
to  the  times  and  positions  of  interest.  For  many  atmospheric  and  ionospheric 
applications  it  is  not  merely  sufficient  to  measure  the  total  electron  fluxes, 
but  the  energy  spectra  and  pitch  angle  distributions  must  be  measured  as  well. 

The  electrons  at  satellite  altitudes  that  are  of  interest  here  are 
those  having  pitch  angles  within  the  loss  cone.  At  middle  latitudes  a large 
fraction  of  the  electrons  observed  at  altitudes  of  a few  hundred  kilometers 
are  generally  locally  trapped  and  therefore  of  no  direct  concern.  But,  the 
distortion  of  the  earth's  magnetic  field  is  such  that  the  altitude  of  the 
mirroring  electrons  changes  significantly  as  the  electrons  drift  eastward 
in  longitude.  The  mirroring  points  dip  to  their  lowest  altitudes  in  the 
vicinity  of  the  South  Atlantic  Anomaly.  Those  electrons  having  pitch  angles 
such  that  the  mirroring  altitude  dips  below  ~ 100  km  are  lost  into  the 
atmosphere.  Since  trapped  electrons  are  continually  being  perturbed,  some 
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fraction  of  those  trapped  at  higher  altitudes  are  continually  being  scattered 
into  lower  mirroring  altitudes  (imhof  and  Smith,  1965) . This  buildup  with 
longitude  drift  of  quasi-trapped  electrons  and  their  subsequent  precipitation 
generally  in  the  longitude  range  from  68°E  to  90%  is  often  called  the 
"windshield  wiper"  effect.  Although  many  of  these  electrbns  are  lost  into 
the  atmosphere  on  the  western  edge  of  the  Anomaly  the  losses  do  occur  on  a 
worldwide  basis.  In  addition  to  the  precipitation  associated  with  a gradual 
change  in  pitch  angle,  electrons  may  be  scattered  more  dramatically  and  be 
lost  into  the  atmosphere  on  a single  bounce.  For  any  of  these  classes  of  pre 
cipitation,  to  obtain  the  energy  deposited  into  the  atmosphere  as  a function 
of  altitude  it  is  necessary  to  measure  both  the  pitch  angle  distributions 
and  the  energy  spectra.  From  such  measurements  the  altitude  profile  of  en- 
hanced ionization  can  then  be  derived. 

As  a result  of  the  Office  of  Naval  Research/Defense  Nuclear  Agency 
sponsored  program  for  deriving  D-region  effective  recombination  rate  coef- 
ficients (Contract  N0001U-70C-0203) , valuable  data  on  the  general  subject  of 
electron  precipitation  at  middle  latitudes  emerged.  These  findings  (Larsen 
et  al.,  1976a)  show  that  energetic  electrons  may  precipitate  into  the  lower 
ionosphere  for  up  to  10  days  following  certain  magnetic  storms  and  thus  very 
likely  may  be  the  most  important  cause  of  the  so-called  "magnetic  storm 
after-effects"  (or  "post-storm  effects")  which  occur  at  middle  latitudes. 

The  electron  precipitation  may  be  strongly  dependent  upon  latitude  as 
shown  in  Figure  2-2  (Larsen  et  al. , 1976b).  For  the  two  ionospheric  stations 
Ottawa  and  University  Park,  Pa.  the  D-region  ionization  may  on  occasions 
very  different,  see  Figures  2-3a  and  2-3h,  albeit  the  close  proximity  in 
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Figure  2-2.  Integral  counting  rates  for  energies  > 130  keV  as  a function  of 
latitude  at  9O0  (HEES  detector)  and  25°  (REES  detector)  pitch 
angles  with  the  geomagnetic  field.  The  University  Park  and 
Ottawa  positions  are  marked  by  arrows.  Panel  (a)  is  for  a pass 
on  16  December  1971  and  (b)  for  a pass  on  19  December  1971 
(Larsen  et  al.,  1976b). 
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ION  PAIR  PRODUCTION  (CM'V1) 

Figure  2-3.  Ion  pair  production  rates  as  a function  of  height  for  the 

University  Park  and  Ottawa  positions,  a)  on  16  December  1971 
and  b)  on  19  December  1971. 
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latitude  (45°  and  4l°N,  respectively)  of  the  sites. 

2.3  The  Satellite  Data  for  Coordinated  Cases 

In  the  previous  report  (imhof  et  al.,  1976)  the  LPARL  satellite  data 
were  presented  and  analyzed  for  11  cases  in  conjunction  with  observed  ELF 
signal  strength  anomalies  over  the  path  WTF  to  Connecticut.  The  UT  times 
and  geomagnetic  conditions  for  these  coordinations  are  listed  in  Table  2-3. 


Table  2-3 


List  of  Cases  for  which  Satellite  Data  were  coordinated  with  ELF 
Signal  Strength  Measurements  (WTF  to  Connecticut)  in  the 
Previous  Effort 


Satellite 

Rev. 

Date 

UT 

(Approx.) 

KP 

I972-O76B 

421 

31  October  1972 

2310 

5“ 

1972 -076b 

422 

1 November  1972 

0049 

8" 

1972-076B 

423 

1 November  1972 

0226 

' 8" 

1972-0 76B 

424 

1 November  1972 

0405 

8 

1971-089A 

5446 

1 November  1972 

0405 

8 

1972-0 76B 

425 

1 November  1972 

0544 

8 

1971-089A 

5447 

1 November  1972 

0545 

8 

1972-076B 

2880 

20  April  1973 

0326 

5 

1972-076B 

3242 

15  May  1973 

044l 

5" 

1972-076B 

6941 

26  January  1974 

0422 

4+ 

1972-0 76b 

7779 

26  March  1974 

04ll 

4" 
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Since  the  time  of  acquisition  of  the  last  data  covered  in  the  pre- 
vious report  extensive  coordinations  were  performed  in  July -August  1975  and 
in  March-April  1976.  During  the  latter  coordinated  exercise  with  Dr.  John 
Davis,  of  NRL,  excellent  satellite  coverage  was  obtained  over  all  of  the  ELF 
stations  in  operation  (Maryland,  Greenland,  Norway  and  Italy) . These  data 
were  acquired  from  the  satellite  1972-076B.  A tabulation  of  the  satellite  data 
acquisition  times  is  provided  in  Table  2-4,  along  with  a listing  of  the  ELF 
receiving  stations  from  which  signal  strength  data  have  been  acquired. 

The  satellite  particle  data  from  all  of  the  passes  listed  in  Table 
2-4,  were  surveyed  with  a special  alert  for  any  unusual  behavior  in  the  elec- 
tron precipitation  profiles,  i.e.  high  fluxes  associated  with  relatively  broad 
precipitation  regions.  Special  attention  was  also  devoted  to  passes  occurring 
at  times  when  unusual  behavior  was  observed  in  the  ELF  data.  For  the  more 
important  passes  the  detailed  energy  spectra  were  processed  and  averaged  over 
broad  latitude  intervals,  energy  deposition  profiles  obtained,  and  electron 
density  profiles  derived.  Selected  portions  of  these  data  are  presented  in 
the  various  sections  of  this  report  dealing  with  the  different  coordinations. 
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Section 


THE  ELF  PROPAGATION  DATA  AND  CALCULATION 
OF  PROPAGATION  CHARACTERISTICS 


3.1  The  ELF  Propagation  Measurements 

The  ELF  field  strength  measurements  were  made  using  signals  from  the 
transmitter  at  the  U.  S.  Navy  ELF  Wisconsin  Test  Facility  (WTF) . The  trans- 
mission source,  located  in  the  Cheqeamegon  National  Forest  in  north-central 
Wisconsin,  consists  of  two  22.5  km  north-south  antennas  (one  buried  and  one 
elevated)  and  one  22.5  km  elevated  east-west  antenna.  Each  antenna  is 
grounded  at  both  ends.  The  transmission  station  is  located  at  the  mid-point 
intersection  of  the  two  antennas.  The  antenna  array  pattern  can  be  steered 
to  any  particular  receiving  location. 

A separate  experimental  ELF  recording  program  under  the  auspices  of 
the  Office  of  Naval  Research  has  provided  measured  field  strength  values  from 
a number  of  receiving  sites  using  WTF  transmissions.  In  this  report  emphasis 
will  be  given  to  the  path  from  WTF  to  Tromso,  Norway;  but  transmission  from 
WTF  to  other  stations  at  Maryland,  Connecticut,  Greenland,  Norway  and  Italy 
will  also  be  discussed.  The  ELF  data  for  the  polar  path  to  Norway  were  ob- 
tained by  the  Naval  Research  Laboratory  under  the  direction  of  Dr.  John  Davis. 
Narrow  band  ELF  receivers  with  effective  integration  times  of  ~ 4000  second 
were  employed.  The  data  were  taken  at  42  or  76  Hz  and  represent  the  actual 
measured  field  strengths  for  two  periods  in  1975  and  one  in  1976. 
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3.2  The  ELF  Computer  Program  for  Model  Waveguide  Calculations 

The  waveguide  model  computer  program  used  in  the  present  studies 
was  developed  at  Naval  Electronics  Laboratory  (NELC)  [ now  Naval  Ocean 
Systems  Center  (NOSC) ] and  is  well  documented  (see  e.g.  Sheddy  et  al.,  1968; 
Pappert,  1970;  Pappert  and  Moler,  1974).  Mr.  W.  Moler  of  NOSC  has  kindly  as- 
sisted in  adapting  the  program  to  run  at  the  Lockheed  (IMSC)  computer  facili- 
ties . 

The  program  is  based  upon  the  formalism  developed  by  Budden  (1961) 
and  Budden  and  Daniell  (1965)  for  electromagnetic  wave  propagation  in  the 
terrestrial  waveguide.  Special  emphasis  has  been  given  to  its  applicability 
in  the  extremely  low  frequency  (ELF)  range. 

The  program  allows  for  a terrestrial  waveguide  with  the  following 
characteristics : 

• Arbitrary  electron  and  ion  density  distributions  with  height. 

• Variable  electron  and  ion  collision  frequency  height  profiles. 

• The  lower  boundary  of  the  waveguide  is  assumed  to  be  a smooth 
homogeneous  earth  with  arbitrary  conductivity  and  permittivity. 

• Allowance  is  made  for  the  anisotropy  of  the  ionosphere;  arbi- 
trary values  for  the  strength  and  direction  of  the  earth’s 
magnetic  field  may  be  applied. 

• Horizontal  inhomogeneties  (along  the  propagation  path)  may  be 
taken  into  account  by  segmenting  the  path  into  regions  having 
different  values  for  the  waveguide  characteristics  listed  above. 

The  computer  program  is  thus  capable  of  treating  waveguides  that  are 
reasonable  model  representations  of  the  actual  complex  terrestrial  waveguide. 
The  program  does  not  allow  for  variation  of  the  guide  parameters  perpendicular 
to  the  propagation  path.  This  simplification,  acknowledged  by  Pappert  and 
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Moler  (197*0  as  the  major  assumption  in  the  model,  should  be  recognized. 

Another  computer  program  (also  developed  at  NELC  and  supplied  to 
LMSC)  is  used  to  calculate  the  signal  strength  at  the  ELF  receiver  using 
data  provided  by  the  waveguide  program  as  input.  The  program  allows  for  seg- 
mentation of  the  propagation  path. 

lor  ELF  propagation  in  the  earth-ionosphere  waveguide  the  expression 
for  the  magnetic  field  component  at  a distance  C from  a horizontal  dipole 
launching  end-on  is  derived  from  expressions  given  by  Pappeyt  and  Moler  (1974) 


H*  «>  - - 


Ez  <«> 


« cowl  /X?  T(o)  E (o)  * exp 

Is H7/.~  k ex  ec 


where 


E^,  is  the  vertical  electric  wave  field 

S is  the  sine  of  the  complex  angle  of  incidence 

Zq  is  the  free-space  impedance 

<t  is  the  azimuth  angle  in  degrees 

a is  the  earth's  radius  in  meters 

T R 

X , X are  the  excitation  factors  for  a horizontal  dipole  launching  end- 

on  at  the  transmitter  and  receiver,  respectively,  at  ground  level 
T R 

ex(o),  e ^ (o)  are  the  height  gain  factors  at  the  transmitter  and 
receiver,  respectively. 

The  factor  Q = 9-01  x 10-8  ILf3/2  [p,V/m]  where 
I is  the  input  antenna  current  in  amperes 

L is  the  length  of  the  transmitting  antenna  in  meters,  and 

f is  the  frequency  transmitted  in  Hz. 
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For  an  antenna  at  ground  level  it  can  be  shown  that: 


eT(o) 


,/ct  /iiue 

V T < 


where 


is  the  earth 
(°EW  = 3-2  X 


conductivity  under  the  transmitting  antenna 
10-^mho/m  for  the  EW  antenna) 


U) 


and 


is  the  angular  frequency  of  the  radiowaves 
(u>  = 2nf) 


is  the  permittivity  of  free  space. 


In  order  to  model  the  ambient  nighttime  conditions  the  following 
electron  and  ion  density  profiles  were  used: 

0 - 100  km 


100  - 200  km 


GE  Tempo,  RPT66TMP-82 
1 March  1 966  by  W.  Knapp 

Satellite  Environment  Handbook, 
Stanford  University  Press,  1965> 
edited  by  F.  S.  Johnson 


The  profiles,  taken  from  pappert  and  Moler  (197*0,  are  shown  in 
Figure  3-1- 

For  the  disturbed  conditions  the  ambient  electron  density  profile  has 
been  modified  above  30-50  km  in  order  to  provide  a smooth  transition  to  the 
N(h)  profile  deduced  from  the  satellite  data  for  greater  heights.  The  N+ 
profile  has  likewise  been  modified  as  discussed  in  later  sections. 


3-4 

LOCKHEED  PALO  ALTO  RESEARCH  LABORATORY 
1 o < ■ n 1 1 • *1111111  » iftci  tontm,  ihc. 

v ? . *; «? *11  « 


atriKNCli: 

o-ioo«m.g(  timio  m.  mtmi-ij. 

IMAICH  1966,  »y  W K N All 

(OO-XOIM,  SATflllTI  INVUONMfNT 

HANMOOI,  STANJOIO  UNIV. 

mu.  i»*i  icxhoiy 

I V JOHNSON. 


IMSC/D560323 


In  this  study  the  nighttime  collision  frequencies  used  by  Pappert  and 
Moler  (197*0  have  been  applied.  The  values  for  the  electron  and  ion  collision 
frequencies  used  are  listed  in  Table  3-1* 


J 


J 

J 


For  propagation  from  WTF  to  Tromso,  Norway  with  path  length  approx 


imately  6 Megameters  (Mm) , we  have  computed  the  radiowave  fields  as  a function 
of  distance  due  to  the  direct  wave  only  and  rejected  the  contribution  from 
the  wave  that  propagates  around  the  world  through  the  antipode  to  Tromso. 

The  difference  in  propagation  length  between  the  long  and  short  path  is  ap- 
proximately 34  Mm.  Therefore,  only  small  inaccuracies  are  introduced  by 
neglecting  the  long  path  component  to  Tromso. 

The  propagation  characteristics  of  ELF  waves  are  dependent  not  only 
upon  the  ionospheric  conditions  along  the  path,  but  also  upon  changes  in  the 
earth's  conductivity  and  permittivity.  The  path  to  Norway  exhibits  large 
changes  in  these  ground  parameters  with  segments  of  sea  water  as  well  as  a 
<*>  1.1  Mm  segment  over  the  Greenland  ice  cap,  as  shown  in  the  map  in  Figure 
3-2.  In  order  to  take  these  changes  into  account,  we  have  divided  the  path 
into  8 segments  which  we  feel  is  a realistic,  yet  practicable  number  of  seg- 
ments to  deal  with.  This  segmentation  is  also  indicated  in  Figure  3-2,  and 
Table  3-2  shows  the  various  characteristic  parameters  of  the  segments.  The 
earth's  magnetic  field  strengths  are  given  for  ground  level  and  these  values 
together  with  the  magnetic  co-dip  angles  and  propagation  azimuths  (with 
respect  to  magnetic  north)  represent  values  for  approximately  the  mid-points 
of  each  segment.  The  corresponding  values  for  the  earth's  conductivity 
(o)  and  relative  permittivity  (e/eQ)  are  supplied  from  Mr.  Moler  of  NOSC. 

For  some  segments,  especially  segment  number  4,  referred  to  as  North  Canada, 
the  a and  e/eQ  values  are  not  constant.  Representative  values  have  been 
selected  for  these  segments.  In  the  ELF  computations  the  changes  in  the  a 
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Segmented  Propagation  Path  WTF 
To  Tromso 


IMSC/D560323 


Field  strength  calculations  have  been  made  for  the  path  to  Tromso 

and  e/e  values  were  modeled  to  take  place  over  a 100  km  path  length  between 
o 

different  segments. 

During  the  observational  periods  of  interest  here  both  the  north- 
south  (NS)  and  east-west  (EW)  antennas  were  activated  at  WTF.  Relative 
phasing  between  the  two  antennas  were  0°,  and  transmitted  frequency  was  76 
Hz  except  during  March/April  1975  when  42  Hz  signals  were  transmitted. 

The  antenna  configuration  for  propagation  to  Tromso,  Norway  is  shown 
in  Figure  3-3,  which  also  illustrates  the  fact  that  the  propagation  direction 
to  Tromso  is  perpendicular  to  the  electrical  axis  for  the  EW  antenna.  At 
Tromso,  therefore,  only  signals  from  the  NS  antenna  are  received  and  the 
signal  strength  should  therefore  be  independent  of  the  NS-EW  antenna  current 
phasing. 

Measurements  at  several  recording  sites  indicate,  however,  that  the 
effective  dipole  moment  of  the  antenna  changes  with  different  phasing  angles 
(Bannister,  private  communication,  1976)  . In  Table  3-3,  for  selected  phase 
values  the  corrections  are  listed  (in  dB)  @ 75  Hz  that  should  be  added  to 
field  strength  predictions  when  the  normalized  antenna  dipole  moment  of 
6.75  x 10  Am  (antenna  current  300A,  antenna  length  22.5  km)  and  = 

3.2  x 10~Snho/m  are  used.  In  our  presentation  both  uncorrected  values  take 
into  account  the  10°  offset  between  the  propagation  direction  to  Tromso  and 
the  direction  of  the  NS  antenna  electrical  axis  (-0.1  dB  in  signal  strength). 
All  field  strengths  throughout  this  report,  except  those  for  the  hypothetical 
Es  layer  are  based  on  a = 3.0  x 10'\nho/m  and  would  be  0.3  dB  lower  for 
a = 3.2k  x 10  ^mho/m. 

No  data  are  presently-  available  as  to  the  corrections  that  should  be 
applied  to  predictions  made  for  42  Hz. 
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ELECTRICAL  AXIS  OF 


. WTF  antenna  configuration  for  propagation  to  Norway. 
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Field  strength  calculations  have  been  made  for  the  path  to  Trctnso 
using  the  set  of  ambient  nighttime  electron  and  ion  density  profiles  as  well 
as  electron-neutral  and  ion-neutral  collision  frequencies  described  in  the 
previous  section.  The  derived  propagation  constants  @ 42  and  75  Hz  are  shown 
in  Tables  3-4  and  3-5,  respectively.  The  values  are  listed  for  each  segment. 
It  is  noted  that  the  Greenland  segment,  as  expected  shows  a higher  at- 
tenuation rate  and  wave  excitation  factor  than  the  other  segments. 

Using  the  propagation  constants  presented  in  Tables  3-4  and  3-5, 
the  ambient  nighttime  field  strength  for  the  path  WTF  to  Tromso  has  been  cal- 
culated. The  results  are: 

42  Hz  75  Hz 

-159.3  dB  -155-6  dB 

Signal  strengths  are  with  respect  to  1 ampere /meter . (A/m) 

No  corrections  for  antenna  phasing  pattern  have  been  applied  to 
these  values.  With  correction  for  0°  phasing  the  corrected  value  at  75  Hz 
becomes:  -154.7  dB  with  respect  to  1 A/m 

3-4  Calculations  for  the  Propagation  Path  WTF  to  Italy,  Maryland, 

Connecticut  and  Greenland 

With  similar  procedures  the  ambient  nighttime  field  strengths  have 
been  calculated  for  the  ELF  receiving  stations  at  Italy,  Maryland,  Connecticut 
and  Greenland.  For  each  of  these  paths  the  segmentation  and  characteristic 
parameters  of  the  segments  are  provided  in  Tables  3-6,  3-7,  3-8  and  3-9. 

Using  the  propagation  constants  presented  in  the  aforementioned 
tables  the  ambient  nighttime  field  strengths  for  the  paths  from  WTF  to  the 
various  stations  have  been  calculated  and  are  summarized  in  Table  3-10. 
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Propagation  Characteristics  @ 42  Hz  for  the  Segmented  Path  for  WTF  to  Tromso,  Norway 


Propagation  Characteristics  @ 75  Hz  for  Segmented  Path  WTF  to  Tromso,  Norway 
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Section  4 

THE  SATELLITE  AND  MULTI-STATION 
COORDINATED  CASES  IN  MARCH-APRIL  1976 

In  March  and  April  1976  a special  coordinated  exercise  was  performed 
between  the  1.PARL  satellite  payload  and  the  ELF  receiving  stations  at  Mary- 
land, Greenland,  Norway  and  Italy.  For  this  coordination  the  payload  was 
operated  at  all  practical  times  when  the  satellite  passed  through  regions 
near  the  transmission  path  and  the  transmitter  and  receivers  were  in  opera- 
tion. The  coordination  with  the  ELF  receivers  was  performed  through  Dr. 
John  Davis  of  the  Naval  Research  Laboratory. 

Precipitating  electron  data  from  the  satellite  passes  during  the 
March-April  1976  coordinations  were  surveyed,  plots  were  made  of  the  pre- 
cipitating fluxes,  and  the  energy  spectra  from  the  more  interesting  passes 
were  obtained  and  subsequently  processed  through  the  LMSC  energy  deposition 
program,  AURORA.  Maximum  precipitating  electron  fluxes  measured  on  each 
satellite  pass  are  plotted  in  Figures  4-1  and  4-2  as  a function  of  time 
superposed  on  the  ELF  field  strengths  measured  at  the  various  receiving 
stations.  The  ELF  transmission  data  were  taken  predominately  on  the  night- 
side.  For  reference,  the  time  of  local  midnight  at  each  station  is  marked 
with  an  arrow.  The  fluxes  of  electrons  > 160  keV  are  shown  in  the  former 
figure  whereas  the  lower  energy  ( > 1.03  keV)  fluxes  are  plotted  in  the 
latter  graph.  To  provide  an  indication  of  the  geomagnetic  activity  plots 
of  the  Dst  index  are  also  included.  Certain  obvious  features  in  the  plots 
are  to  be  noted.  The  various  receiving  stations  display  a similar  nightly 
average  pattern,  i.e.  a general  trend  toward  increasing  signal  strength  as 

4-1 

LOCKHEED  PALO  ALTO  RESEARCH  LABORATORY 
(•ciniii  a 1 1 1 hi  1 a tract  coaranv.  inc. 


4-2.  Plot  of  the  precipitating  electron  flux  at'  1 keV  as  measured  from  the  satellite 
on  coordinated  passes  and  the  ELF  signal  strengths  as  received  at  the  various 
stations.  For  comparison,  the  Ds-fc  geomagnetic  index  is  also  plotted  to  indicate 
the  general  occurrence  of  magnetic  storms. 
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time  progresses  from  March  23  onward  with  a leveling  off  at  a broad  maximum 
in  early  April  followed  by  a decrease  toward  the  middle  of  April.  On  a broad 
time  scale  the  various  stations  indicate  a consistent  pattern  which  may  be 
associated  with  the  geomagnetic  storms  on  26  March  and  1 April.  On  the  other 
hand,  when  one  examines  the  signal  strength  profiles  on  a shorter  time  scale 
it  is  clear  that  large  inconsistencies  exist  between  the  various  stations. 

All  of  the  measured  variations  may  be  due  to  the  effects  of  energetic  particle 
precipitation  since  the  latter  phenomena  are  known  to  be  highly  variable  in 
space  and  time.  For  these  reasons  the  precipitating  electron  fluxes  and 
energy  spectra  should  be  measured  simultaneously  over  a broad  range  of  longi- 
tudes. This  would  involve  the  use  of  particle  spectrometers  on  multiple 
satellites  or  a single  x-ray  mapping  spectrometer  on  one  satellite  (imhof 
et  al.,  197^A) . The  latter  technique  has  many  attractive  features  for  the 
present  purposes,  but  such  data  have  yet  to  be  taken  with  adequate  spatial 
resolution  over  the  energy  range  of  interest.  In  view  of  the  complex  nature 
of  the  data  and  the  present  lack  of  detailed  electron  precipitation  maps 
covering  the  broad  range  of  longitudes  involved,  the  principle  efforts  in 
the  present  study  have  been  primarily  concerned  with  nightly  averaged 
signal  strengths. 

Before  considering  the  energetic  particle  data  let  us  first  investigate 
possible  correlations  between  the  ELF  signal  strengths  and  the  geomagnetic 
activity  K index  as  recorded  at  a magnetic  observatory  near  the  ELF  receiving 
station.  Such  a comparison  with  the  K indices  is  performed  for  general 
scoping  purposes  with  the  full  realization  that  the  precipitating  electron 
fluxes  have  been  demonstrated  not  to  be  correlated  in  detail  with  geomagnetic 
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activity.  In  Figures  4-3  and  4-4,  for  stations  at  Tromso  and  Maryland 
respectively,  these  two  quantities  are  plotted  with  overlapping  scales 
against  time  during  the  period  March-April  1976.  The  appropriate  data 
have  been  averaged  over  the  indicated  time  intervals  to  minimize  fine 
time  scale  variations  consistent  with  restriction  to  nightime  measurements. 

At  each  of  these  stations  the  variations  of  the  two  quantities  display 
many  similarities,  suggesting  some  degree  of  correlation.  The  indicated 
similarities  are  in  the  direction  of  acquiring  higher  signal  strengths  when 
more  geomagnetic  activity  occurs  and  hence  when  the  fluxes  of  precipitating 
particles  are  expected  to  be  greater.  The  reader  should  be  cautioned, 
however,  that  despite  the  consistency  displayed  in  Figures  4-3  and  4-4, 
the  agreements  are  less  convincing  when  similar  plots  are  made  for  other 
time  periods.  This  is  illustrated  in  Figures  4-5  through  4-6  for  the  time 
periods  March-April  1975.  Some  features  appear  with  a similar  behavior 
in  both  the  signal  strengths  and  the  K index  over  the  region  of  the  receiver. 
However,  there  are  also  many  examples  of  opposite  behavior.  For  example, 
the  ELF  signal  strengths  at  Maryland  were  below  the  threshold  of  detection 
at  times  in  March  1975  when  intense  relativistic  electron  precipitation 
occurred.  Similar  plots  of  the  coordinated  data  taken  in  July-August  1975 
are  not  presented  here  because  of  the  daylight  nature  of  many  of  the  ob- 
servations. 

A correlative  analysis  was  performed  on  the  basis  of  the  particle  data 
survey  and  the  ELF  data  provided  by  Dr.  Davis.  For  Figure  4-7  the  measured 
ELF  signal  strengths  are  plotted  against  the  precipitating  electron  fluxes 
for  the  Maryland,  Greenland,  and  Norway  stations  in  March-April  1976.  Thf3 
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Figure  4-7.  Measured  ELF  signal  strengths  at  76  Hz  as  received  at  Maryland 
Greenland  and  Norway  stations  plotted  as  a function  of  the 
integral  fluxes  of  precipitating  electrons  above  energies  of 
160  keV  and  the  differential  precipitating  electron  fluxes  at 
1.03  keV  + 0.45  keV. 

4-10 

LOCKHEED  PALO  ALTO  RESEARCH  LABORATORY 

10CINIIS  *1111111  t l»*C«  *•■»*»».  <NC. 


IMS  C/D  560323 


data  points  clearly  display  a largo  amount  of  scatter,  masking  any  general 
trends.  The  measurements  may  suggest  some  tendency  toward  higher  signal 
strengths  at  greater  electron  fluxes,  but  this  is  not  clearly  established 
with  the  present  data  base.  In  some  earlier  Connecticut  data  presented  in 
Figure  4-8  there  appears  to  be  a general  decrease  in  signal  strength  with 
increasing  flux  of  the  precipitating  electrons  above  160  keV.  These  trends 
are  not  suggested  in  the  fluxes  of  1 keV  electrons  which  deposit  their  energy 
in  the  upper  portions  of  the  E-region.  From  these  figures  one  is  left  with 
the  impression  that  energetic  ( > l60  keV)  particle  precipitation  may  effect 
the  ELF  signal  strengths,  but  that  either  attenuations  or  amplifications 
may  occur  as  a result  of  increase  precipitating  flux.  Clearly,  a more 
detailed  analysis  of  the  data,  making  use  of  the  intensity  profiles  and 
energy  spectra  is  needed  to  be  performed  in  conjunction  with  calculations 
of  the  expected  signal  strengths.  In  line  with  these  objectives  the 
satellite  data  have  been  analyzed  in  detail,  as  discussed  below. 

For  several  coordinated  satellite  passes  during  the  period  March-April 
1976  the  measured  fluxes  of  precipitating  electrons  > 160  keV  are  plotted  in 
Figure  4-9  as  a function  of  invariant  latitude.  The  latter  parameter  is  a 
geomagnetic  coordinate  appropriate  for  study  of  energetic  particle  precipi- 
tation. For  reference  the  invariant  latitude  of  the  WTF  transmitter  is  59° 
and  the  receiving  stations  at  Connecticut,  Maryland,  Greenland  and  Norway 
are  57°>  54°,  86°  and  66°,  respectively.  Since  the  waveguide  program  pro- 
vides for  segmenting  the  propagation  path  into  different  regions  but  does 
not  allow  for  variation  of  the  guide  parameters  perpendicular  to  the  propa- 
gation path  some  type  of  latitude  averaging  must  be  invokved,  and  therefore 

4-11 


LOCKHEEO  PALO  ALTO  RESEARCH  LABORATORY 

A""".*  . -L’.ML V/  . J“-: 


ocr.  1972 -MAY  1973 


imm 


IMSC/D560323 


sS  V. 


+ 1 o 

^ w 

_ I 

© 5 

^ a 


» X 

s3 


(DOC 

go 

hO  C 
O C -P 
O *H  O 
-P  o 
(D  CD  H 
X -P  o 
-P  «H 
Q4  bD 
-P  C 

CO  O *H 
O -P 
T3  (0 

(D  CD  -P 

> -rH 

•H  Qj 

o O *H 
O O 

o C/3  0 
0 fn 


w 

"o  J 

S w 


CD  C 
vD  ^ a) 
t'-  bD 
0)  V 
P P 
CD  C <*h 


1 — I "-H  £> 

CO  -P  V <D 


5 'o 

x W 

S T 
s? 
a a 

h a 


0)  T3  <U  4) 
l(  OJ  P JC 
3 -P  bD 
«1  p ti  (D 

CO  O <D  O 
flj  H c • 
2:  Q.  <U  f-l 


(w/v  i ox  nan  ap)  hxonshxs  ivnoio  h 9 l 


* 


800  UT  26  MARCH  1976 


M3C/D560323 


I 

L 


I 


r* 


II 

u 


-L-L 


till  l I I I I I I 

"t.  N. 


<M 

o 


o 

a 

•H 

-P 

•r( 

o 

a; 

co 

•H 

<4H 

o 

Jh 

CO 

<v 

X 

as 

0) 

T3 
T3  <U 
I)  -P 
P>  « 
-P  C 
o -rj 

r-l  T 3 
04  ^ 

O 

> O 
0)  o 

■id 


Q) 

s 

X 

-p 

rH 

bC 

A 

C 

•H 

CO 

H 

c 

3 

o 

TJ 

-p 

CO 

o 

<D 

<D 

CO 

•H 

CO 

CD 

CO 

04 

W 

c 

T3 

•H 

CD 

-P 

■P 

CO 

O 

■P 

CD 

•H 

rH 

04 

<b 

•H 

CO 

O 

d) 

Sh 

O 

04  <P 

Cm 

CD 

O 

T*  • 

3 VO 

CO 

-p  h- 

a> 

•H  ON 

X 

-P  rH 

3 

CO 

rH 

»H  rH 

Qh 

H 

-P  u 

T3 

a 04 

CD 

CO  < 

U 

•H  1 

3 

CO 

CO  o 

CO 

> ^ 

CD 

C CO 

s 

2 

• 

«; 

CT\ 

1 

-J- 

<d 

3 

t>D 

•H 

f*4 

(,_K-  WO)  AM  091  < xnn  0NI1V1MO3W  UN 


4-13 


NET  PRECIPITATING  FLUX  > 160  keV  (CM-2 -SEC*1 -SR 


LMSC/D560323 


0110  6 APRIL  1976 


Figure  4-  9C  . Measured  fluxes  of  precipitating  electrons  > 160  keV  plotted 

as  a function  of  invariant  latitude  for  a selected  pass  during 
the  coordinated  exercise  in  March-April  1976. 
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several  empirical  approaches  were  taken.  These  included  a linear  averaging 
of  the  energy  deposition  profiles  over  the  complete  magnetic  latitude 
interval  between  the  transmitter  and  the  receiving  station  as  well  as  other 
procedures  such  as  selecting  the  highest  fluxes  within  an  interval.  For 
many  of  the  analyses  presented  here  maximum  fluxes  over  a latitude  interval 
of  at  least  a few  tenths  of  a degree  were  used. 

Examples  of  the  more  intense  energy  spectra  of  precipitating  electrons 
measured  during  the  March-April  1976  coordinations  are  shown  in  Figure  4-10. 

These  measurements  were  all  performed  near  local  midnight,  and  in  each  case 
all  or  virtually  all  of  the  propagation  path  of  interest  was  in  darkness. 

Overall  the  precipitating  electron  intensities  are  quite  high  but  somewhat 
lower  than  many  recorded  during  the  31  October-1  November  1972  magnetic 
storm  event  when  high  fluxes  of  precipitating  electrons  were  observed  and 
when  pronounced  ELF  transmission  anomalies  were  recorded.  The  spectra 
shown  indicate  significant  intensities  at  energies  > 1 MeV,  but  spectra 
with  much  greater  fluxes  of  relativistic  electrons  have  been  observed  and 
are  presented  in  Section  5. 

For  each  of  the  electron  energy  spectra  processed  the  ion-electron 
production  rates  in  the  atmosphere  were  calculated  as  a function  of  altitude, 
using  the  computer  program  AURORA  (Vvalt  et  al.,  1968) . The  method  consists 
of  solving  the  time-independent  Fokker-Planck  diffusion  equation  of  elec- 
trons by  use  of  finite  difference  techniques.  The  mirroring  effects  of  the 
converging  magnetic  field  are  included.  The  program  accepts  an  input 
electron  differential  energy  spectrum  of  arbitrary  shape  along  with  any 
pitch  angle  distribution.  Energy-dependent  pitch  angle  distributions  may 
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also  be  applied.  In  such  a case  the  program  uses  the  input  pitch  angle 
distributions  prescribed  for  selected  energies  to  obtain  distributions  at 
other  energies  through  interpolation.  In  the  calculations  the  electrons 
are  allowed  to  interact  with  the  atmosphere  below  300  km.  The  program 
outputs  include  the  electron  pitch  angle  distribution  between  0°  and  l80°, 
thus  including  backscattered  and  reflected  electrons,  at  pre-selected 
energies  and  altitudes  as  well  as  the  total  energy  deposition  within  speci- 
fied altitude  regions.  Examples  of  ion-electron  production  rate  profiles 
are  shown  in  Figure  4-11. 

The  ion-pair  production  rates  shown  in  Figure  4-11  are  in  general  quite 

large  compared  to  the  natural  ionization  rates  from  cosmic  rays.  For  example, 

at  70  km  the  cosmic  ray  production  rates  are  less  than  0.1  ion  pair/cm^-sec. 

The  ion  pair  production  rates  encountered  in  the  events  of  present  interest 

are  to  be  compared  with  those  obtained  during  the  very  intense  solar  particle 

events  in  August  1972  (Reagan  and  Watt,  1976).  During  the  peak  of  the  August 

1972  event  the  production  rates  over  the  Chatanika,  Alaska  radar  site  were 

less  than  l(r  ion  pair/cm  -sec  at  all  altitudes  below  90  km.  Of  course,  due 

to  the  high  energy  protons  in  the  PCA  event  the  ion  pair  production  rates 

were  much  higher  at  the  lower  altitudes.  For  example,  at  40  km  over 

Chatanika  at  1508  UT  on  4 August  1972  the  production  rate  was  about  3 x 10^ 

3 

ion-pair/cm  -sec.  However,  at  an  altitude  of  90  km  the  production  rate  was 
3 3 

only  about  3 x Hrion-pairs/cm  -sec  which  is  to  be  compared  with  the  rates 
in  Figure  4-11  that  are  higher  by  one  to  two  orders  of  magnitude.  This 
difference  largely  reflects  the  higher  fluxes  of  electrons  in  the  tens  of 
keV  energy  range  during  the  events  covered  in  Figure  4-11.  Likewise,  during 
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r production  rates  obtained  from  the  measurements  of  precipitating  electrons 
selected  satellite  passes  in  March-April  1976. 


WATT,  NEWKIRK,  AND  SHELLEY  (1974) 


Figure  4-12,  Effective  electron  loss  rates  plotted  as  a function  of  altitude.  D-region 
data  are  shown  from  the  polar  cap  absorption  events  on  19  November  1968 
and  3 November  1969  (Ulwick)  and  from  the  PCA  event  in  August  1972  (Reagan) . 
At  higher  altitudes  data  from  Watt  et  al  are  included.  The  selected 
best-fit  curve  is  also  shown. 


some  other  selected  solar  particle  events  such  as  April  13,  1969  Potemra 

o 4 o 

and  Zmuda  (1972)  calculated  production  rates  of  10  -*  10  ion  pairs/cm-sec 

in  the  70  to  9°  ^ altitude  range. 

In  order  to  obtain  electron  densities  from  the  ion  pair  production 
rates  it  is  necessary  to  invoke  values  for  the  effective  nighttime  loss 


rate  coefficients.  Various  measurements  of  the  loss  rate  coefficients 
have  been  made  in  the  past,  but  at  present  there  is  not  a general  concur- 
rence on  the  values.  The  effective  loss  rates  are  known  to  vary  with  the 
production  rate  at  night  in  contrast  to  the  daytime  loss  rates  (Swider  and 
Keneshea,  1972).  At  the  times  of  the  ELF  transmission  anomalies  of  interest 
here  the  ion  pair  production  rates  are  generally  comparable  to  those  at  the 
times  of  PCA  events.  The  loss  rates  obtained  from  PCA  observations  should 
therefore  be  applicable.  These  are  shown  in  Figure  4-x2.  The  loss  rate 
coefficients  are  for  nighttime  conditions  ( x = 98°)  during  the  PCA  event 
in  August  1972  (Reagan  and  Watt  1976) . At  higher  altitudes,  data  from 
Watt  et  al.,  1974  are  also  included.  Based  on  the  published  values  a best- 
fit  curve  was  selected,  as  shown. 

For  analysis  of  these  data  we  have  assumed  the  existence  of  a steady- 
state  condition  in  the  D-region.  Such  an  assumption  simplifies  the  analysis, 
but  requires  that  the  energy  input  not  change  substantially  during  typical 


11 


t* 


recombination  time  periods  which  are  of  the  order  of  tens  of  seconds.  From 
the  data  acquired  on  a single  satellite  pass  it  is  difficult  to  separate 
temporal  and  spatial  events.  The  data  acquired  on  subsequent  passes  (about 
100  minutes  later)  typically  indicate  substantially  different  electron  pre- 
cipitation profiles  but  these  differences  may  largely  reflect  longitude 
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variations.  From  a large  body  of  satellite  data  acquired  by  the  LPARL  group 
on  the  rates  of  precipitation  of  energetic  electrons  over  a large  area  using 

the  atmospheric  bremsstrahlung  technique  (imhof  et  al.,  1974a)  it  has  been 
found  that  the  total  electron  influxes  typically  remain  rather  constant 
during  the  time  ( ~ 5 minutes)  for  traversal  of  the  satellite  over  a given 
region.  The  conclusion  drawn  from  the  satellite  bremsstrahlung  measurements 
should  be  applicable  to  the  present  analyses  where  we  are  primarily  concerned 
with  the  ionization  produced  over  a rather  large  area.  For  a steady-state 
condition  in  the  D-region  the  relation  between  production  and  loss  of  elec- 
irons  may  be  written: 

Q = i|(N2 

where  Q = total  ion  pair  production  rate,  \|i  = effective  electron  loss  rate 
and  N = electron  density.  Based  on  this  relation  and  the  values  of  Q and 
ill  discussed  above  the  electron  densities  were  obtained  as  a function  of 
altitude.  Examples  of  electron  density  profiles  for  the  March-April  1976 
coordinated  exercise  are  shown  in  Figure  4-13. 

At  altitudes  above  ~ 80  km  the  positive  ion  density  profiles  were  taken 
to  be  equal  to  the  electron  densities.  This  altitude  limit  was  somewhat 
- itrarily  selected,  and  may  indeed  be  as  high  as  85-90  km,  but  the  computer 
runs  indicated  a lack  of  sensitivity  to  the  N+  densities  in  this  altitude 
region.  At  altitudes  below  ~ 40  km  ambient  densities  for  the  N+  profiles 
were  used.  In  between  these  altitude  limits  a smooth  curve  interpolation 
was  assumed.  For  the  electron  precipitation  at  0602  HR  UT  on  26  March  1976 
the  electron  and  positive  ion  density  profiles  ar’  shown  in  Figure  4-l4. 


I 


Examples  of  electron  density  profiles  obtained  from  the  March-April  1976  coordinated 
satellite  passes. 
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The  electron  density  profiles  obtained  by  the  above  procedures  were 
then  incorporated  into  the  ELF  propagation  computer  program.  Since  at  any 
given  time  the  satellite  measurements  of  the  precipitating  electrons  were 
performed  at  only  one  longitude  crossing  it  was  necessary  to  make  certain 
assumptions  regarding  the  spatial  extent  of  the  enhanced  ionization  densities. 
The  geometry  of  the  transmitter,  the  satellite  crossing,  the  electron  precipi- 
tation oval  and  the  receiver  are  illustrated  schematically  in  Figure  2-1.  The 
guidelines  used  for  calculation  purposes  were  to  assume  that  the  maximum  pre- 
cipitating electron  inputs  averaged  over  a 4°  wide  latitude  band  applied  over 
the  full  longitude  interval  of  concern  and  over  a latitude  dimension  consistent 
with  the  segmentation  used  in  the  ELF  computer  program.  Even  if  available, 
more  detailed  electron  density  profiles  could  not  be  used  in  the  present 
waveguide  mode  program  due  to  the  important  limitations  of  that  program  in 
not  being  able  to  treat  variations  of  the  guide  parameters  perpendicular  to 
the  propagation  path.  Clearly,  for  future  more  quantitative  studies  this 
major  simplification  in  the  program  should  be  rectified  and  more  detail 
spatial  mappings  of  the  precipitating  electrons  obtained.  The  effect  of  the 
present  simplified  analysis  is  perhaps  to  set  an  upper  limit  to  the  signal 
attenuations,  although  for  any  given  case  the  effect  on  the  ELF  signals  might 
be  either  underestimated  or  overestimated  due  to  the  possible  variations 
with  longitude  in  the  precipitating  electron  fluxes. 

Several  runs  of  the  ELF  waveguide  mode  computer  program  were  made  to 
study  the  effects  on  the  field  strength  of  electron  density  profiles  measured 
during  the  March-April  1976  coorindated  exercise.  The  results  of  these  runs 
for  a representative  electron  spectrum  (0602  UT  on  26  March  1976)  are  listed 
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Section  5 

THE  EFFECTS  OF  RELATIVISTIC  ELECTRON 
PRECIPITATION  EVENTS 

On  several  satellite  passes  during  the  coordinated  exercises  in 
March-April  1975  and  July  1975  intense  relativistic  electron  precipitation 
events  (REP's)  with  very  high  fluxes  of  energetic  electrons  were  observed. 

An  example  of  an  REP  spectrum  measured  on  11  March  1975  is  shown  in  Figure 
5-1  for  comparison  with  other  spectra.  During  the  REP  event  the  intensities 
at  energies  of  several  hundred  keV  and  above  are  much  higher  than  observed 
during  the  31  October-1  November  1972  solar  particle  event  when  high  fluxes 
of  precipitating  electrons  were  measured.  In  the  REP  one  should  note  the 
broad  peaking  near  1 MeV  which  leads  to  a high  energy  deposition  rate  and 
hence  a high  electron  concentration  at  low  altitudes. 

Other  examples  of  REP  spectra  measured  during  the  coordinated  exercises 
are  shown  in  Figure  5-2.  Spectra  of  this  type,  with  broad  peaks  in  the  several 
hundred  keV  to  few  MeV  region,  have  been  observed  many  times  at  high  lati- 
tudes near  the  radiation  belt  trapping  boundary.  These  precipitationg  spectra 
cause  high  ion  pair  production  rates  at  low  altitudes,  as  shown  in  Figure 
5-3 9 but  their  effects  on  ELF  transmission  may  be  somewhat  reduced  by  the 
rather  narrow  latitude  interval  over  which  they  are  effective.  The  REP 
spectra  considered  here  were  all  measured  near  local  midnight.  So  to  obtain 
the  electron  densities  shown  in  Figure  5-4,  use  was  made  of  the  nighttime 
effective  loss  rates  discussed  in  Section  4. 

An  evaluation  of  the  effects  of  the  measured  REP  events  on  the  ELF 
transmission  signal  strengths  has  been  made  with  use  of  the  derived  electron 
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Examples  of  ion-pair  production  rates  from  REP  events 


Examples  of  enhanced  electron  density  profiles  resulting  from  BEP  events 
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density  profiles  incorporated  into  the  NOSC  waveguide -mode  computer  program. 

For  the  bulk  of  the  computer  runs  the  N+  profiles  at  altitudes  below  4o 

km  were  taken  to  be  the  same  as  the  ambient  profiles  discussed  in  Section  3-2. 

As  discussed  in  Section  4 at  altitudes  Ja  80  km  the  N+  densities  were  set 

equal  to  the  electron  densities.  For  altitudes  between  these  two  limits  the 

N+  densities  were  estimated  using  the  ion  pair  production  rates  for  an  ion-ion 

-73-1  + 

recombination  rate  coefficient  of  10  cm  -sec  . Specifically,  the  N den- 
sities were  computed  from  the  relation  N+  = \J q/o\  , where  a ^ is  the  ion-ion 
recombination  rate  coefficient.  This  relation  follows  from  the  following 
considerations: 

The  differential  equation  for  rate  of  positive  ion  production  is: 
dH+  = q - o?i  N+N_  - cDNeN+ 


For  steady  state  conditions  dN  = 0,  therefore 

T 


q = c*.  + VeN+ 


= WNe  + oDNeN+ 


= NeN+(*D  + X^) 


where  X.  = N /N 


or  q ^e  (orD  + Xo\)  N+t’ 

N+ 

At  altitudes  significantly  below  80  km  where  \ is  large  and 


N » N 


then  q N+  or 
or  N+  = ^q/a? 
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In  order  to  gain  an  understanding  of  the  sensitivity  of  the  calculated  signal 
strengths  to  the  N+  profile,  computer  runs  were  also  made  with  a different 
value  of  the  ion-ion  recombination  rate  coefficient.  For  a value  of 

go  *|^ 

6 x 10  cm  -sec  , shown  in  Figure  5-5,  the  signal  strengths  were  not  sig- 

-7  3 -1 

nificantly  different  from  those  obtained  with  10  cm  -sec  . From  these 
computer  calculations  one  concludes  that  over  the  range  of  "reasonable" 
values  of  the  ion-ion  recombination  rate  coefficient  the  transmitted  signal 
strength  is  not  very  sensitive  to  the  N+  profile. 

For  the  event  on  11  March  1975  the . calculations  were  performed  for  paths 
from  WTF  to  each  of  the  following  receiving  stations;  Norway,  Greenland, 

Italy,  Connecticut  and  Maryland.  Since  for  a given  REP  event  the  measure- 
ments were  limited  to  the  immediate  vicinity  of  the  satellite  path  certain 

assumptions  were  made  as  to  the  world-wide  extent  of  the  electron  precipitation. 

o o 

The  high  electron  fluxes,  measured  over  a latitude  extent  of  ~ 0.5  to  ~ 1.0 
at  one  longitude  crossing  of  the  satellite,  were  assumed  to  extend  over  a 
significant  portion  of  the  transmission  path.  Under  these  assumptions  the 
calculated  signal  strengths  for  a relativistic  electron  precipitation  event, 
on  11  March  1975,  are  provided  in  Table  5-1.  The  calculations  were  performed 
at  the  frequency  42  Hz,  in  effect  at  the  time  of  the  electron  measurements. 

The  computer  program  predicts  signal  enhancements  for  the  stations  at  Mary- 
land, Connecticut,  Italy,  Greenland  and  Norway. 

ELF  transmission  calculations  were  also  performed  for  the  intense  REP 
at  0553  UT  on  21  March  1973.  The  computations  were  applied  to  the  path 
"•«T  to  Trotnso,  Norway  for  a variety  of  assumptions  as  to  the  location 

■ ■ : recipitation  consistent  with  the  outer  belt/auroral  zone  region 
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Figure  5-5.  Tv/o  different  positive  icn  density  profiles  corresponding  to  an  ion-ion  recombination 
rate  coefficient  of  10-7Cm3-sec_]-  (Curve  i)  and  6 x 10“8cm3-sec-:L  (Curve  II)  . 
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illustrated  in  Figure  2-1.  Since  the  precipitation  oval  comes  close  to  the 
vicinity  of  the  transmitter  and  the  receiver  at  Tromso  the  ELF  waveguide 
program  was  run  for  precipitation  near  the  transmitter  above,  near  the 
receiver  alone  and  at  both  locations.  The  results  are  summarized  in  Table 
5-2.  It  is  clear  that  either  signal  increases  or  decreases  may  occur,  de- 
pending upon  the  spatial  extent  and  location  of  the  ionization.  Although 
insufficient  coordinated  data  are  available  for  verifying  the  predicted 
effects  of  these  very  intense  and  hard  relativistic  electron  precipitation 
events  the  calculations  based  on  the  satellite  data  indicate  that  either 
signal  enhancements  or  degradations  may  occur.  The  predicted  effects  of  a 
single  relativistic  electron  precipitation  event  are  not  severe,  but  due  to 
their  frequent  occurrence  they  may  provide  the  opportunity  to  verify  ex- 
perimentally the  predicted  effects  of  more  severe  and  rarely  occurring 
phenomena  such  as  PCA  events.  For  a more  precise  evaluation  of  their  continual 
role  in  regard  to  ELF  propagation,  the  acquisition  and  analysis  of  data  on 
more  coordinated  passes  should  be  performed. 

Electron  precipitation  events  can  produce  some  ionization  at  altitudes 
down  to  ~ 30  km  through  the  production  at  higher  altitudes  of  bremsstrahlung 
x-rays  which  can  penetrate  deeply  into  the  atmosphere.  Examples  of  the 
ionization  produced  by  the  bremsstrahlung  are  shown  in  Figure  5-6,  both  for 
an  REP  event  and  for  softer  electron  precipitation  of  the  type  discussed  in 
Section  4.  As  one  can  see  the  ion  pair  production  rates  from  the  bremsstrahlung 
x-rays  may  exceed  nighttime  ion  production  rates  and  hence  can  modulate  mid 
D-region  ionization. 
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-rays  produced  by  precipitating  electrons. 


For  the  REP  event  on  21  March  1973  ELF  signal  strength  calculations 
were  performed  with  inclusion  of  the  additional  ionization  produced  by  the 
bremsstrahlung  generated  in  the  atmosphere  by  the  precipitating  electrons. 


Effect  of  Bremsstrahlung  Produced  in  Precipitation  Event 
at  0553  ITT  on  21  March  1973  Transmission  from  WTF  to 
Receiving  Station  at  Tromso,  Norway 
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Section  6 

WAVEGUIDE  MODE  CALCULATIONS 
FOR  POLAR  CAP  ABSORPTION  EVENTS 


6.1  Calculations  for  Particle  Inputs  Measured  during  the  4 August  1972 

PCA  Event 

No  ELF  propagation  data  from  WTF  have  been  obtained  at  Tromso,  Norway 
during  polar  cap  absorption  (PCA)  events  (up  to  March  1977) • In  order  to 
estimate  the  effects  of  PCA  events,  calculations  were  performed  with  the 
NOSC  waveguide  computer  program  on  the  basis  of  actual  electron  density 
and  energetic  particle  precipitation  data.  Nighttime  PCA  conditions  for 
this  path  were  simulated  using  nighttime  ionospheric  data  from  the  great 
event  of  4 August  1972.  This  PCA  event  was  the  most  intense  in  regard  to 
proton  precipitation  with  energies  between  20-100  MeV  during  the  present 
solar  cycle.  The  data  were  obtained  during  a coordinated  experiment  with 
LPARL  and  Stanford  Research  Institute  using  the  ground-based  Chatanika  in- 
coherent scatter  facility  and  the  1972-076B  satellite  (Reagan,  1975;  Reagan 
and  Watt,  1976). 

In  Figure  6-1  is  shown  the  electron  density  profile  as  measured  by  the 
Chatanika  radar  between  59  and  100  km  at  4 August  1972  ll44  UT  (Ol44  LT)  for 
a solar  zenith  angle  of  x = 95*3°.  For  the  ELF  calculations  it  is  important 
also  to  know  the  N(h)  densities  below  59  'sm.  For  this  purpose  the  production 
rates  (q)  of  ion  pairs  have  been  combined  with  typical  daytime  and  nighttime 
electron  loss  rates  (y)  to  yield  electron  densities  by  the  following 
relation: 

N =“  q/y 
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It  is  surmised  that  these  values  will  serve  as  upper  and  lower  limits  to 

the  electron  densities  present  at  x = 95.3°« 

The  electron  loss  rates  for  daytime  conditions  (4.0  x 10~^  and 
-4  3-1 

4.0  x 10  cm  -s  at  55  and  50  km,  respectively)  were  taken  from  Reagan 
(1975) ; the  value  at  50  km  being  obtained  by  extrapolation.  The  nighttime 
loss  rates  were  obtained  by  extrapolation  to  lower  altitudes  the  values  de- 
duced by  Sellers  and  Stroscio  (1975)  for  the  November  69  FCA 
-2  -1.  -1 

(1.0  x 10  and  1.0  x 10  cm  ^s  at  55  and  50  km,  respectively). 

The  N(h)  profile  below  59  km  used  in  the  ELF  calculations  is  a "best 

guess"  between  the  established  upper  and  lower  limits  with  preference  given 

to  the  lower,  nighttime  values.  Below  40  km  the  N(h)  profile  joins  smoothly 

with  the  ambient  nighttime  profile.  The  positive  ion  density  for  this  case 

is  derived  as  follows:  above  83  km  N+  is  taken  to  be  equal  to  N the  N+  densi- 

e 

ties  at  80,  70  and  60  have  been  computed  for  the  total  duration  of  this  FCA 
event  by  Gunton  et  al.  (1977),  and  the  values  marked  by  open  squares  repre- 
sent the  computed  densities  at  =“  1144  UT  on  the  day  in  question.  At  35  to 
50  km  the  values  were  approximated  by 

n+  q/o,i 

where  is  the  ion-ion  recombination  rate  coefficient.  This  approximation 

is  valid  at  low  altitudes  where 

o'  N+N  < <*,  N+N~ 
a I 

and  steady  state  conditions  are  assumed  (dN+/dt  = 0) . The  factor  is  the 
electron-ion  recombination  rate  coefficient.  Since  at  these  altitudes  the 

“ <f  ■ -f  O 

electron  densities  N < < N , we  have  N =-  N and  q =»  <y^(N  ) . When  a 

-7  Q -1  + 

constant  value  of  cr  = 1.10  cm  -sec  with  height  is  assumed,  the  N values 
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marked  by  open  circles  are  obtained.  Below  25  km  ambient  densities  for 
N+  are  assumed.  The  N+  values  marked  with  asterisks  at  50,  40  and  35  km 
were  obtained  using  the  "selected’’  ar^(h)  values  listed  by  Cole  and  Pierce 
(1965)  . It  may  be  concluded  that  for  our  computations  the  assumption  of 
a constant  or  with  height  is  sufficiently  accurate. 

The  propagation  characteristics  using  these  N and  N+  profiles  for  the 
various  segments  for  the  WTF  path  to  Tromso,  Norway  are  listed  in  Tables 
6-1  and  6-2  for  4-2  and  75  Hz,  respectively.  It  may  be  noted  that  the  at- 
tenuation rates  are  up  significantly  over  the  ambient  nighttime  values. 
Simultaneously,  the  excitation  factor  has  increased  by  5-6  dB. 

Using  these  propagation  constants  the  calculated  field  strength  at 
Tromso  for  75  Hz  is: 

-I58.9  dB  wrt  1 A/M  (uncorrected  for  antenna  phasing) 
which  is  -3.3  dB  lower  than  calculated  for  ambient  conditions. 

It  may  happen,  however,  as  in  the  present  case,  that  the  region  dis- 
turbed by  the  solar  protons  will  not  cover  the  complete  WTF-Tromso  path. 

At  ~ 1200  UT  on  4 August  1972  the  latitude  cutoff  for  20-100  MeV  proton  pre- 
cipitation was  occurring  at  approximately  L 3-8  or  an  invariant  latitude 
of  59°.  In  fact,  had  this  been  a real  coordinated  case  the  WTF  transmitter 
region  would  have  been  outside  the  PCA  disturbed  area. 

If  we  make  the  assumption  that  the  WTF  segment  is  outside  the  PCA  dis- 
turbance and  has  ambient  nighttime  ionospheric  conditions,  while  the  PCA  is 
occurring  over  the  rest  of  the  path,  the  calculated  field-strength  at 
Tromso  for  75  Hz  is 

-l6l.5  dB  wrt  1 A/M  (uncorrected  for  antenna  phasing) 
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Propagation  characteristics  @ 42Hz  for  the  segmented 


40 

x 

CO 

03 

bD 

0 

T3 

•r-l 

CO 

CO 

CO 

CO 

0 

co 

G 

T3 

C 

G 

cq 

CO 

40 

C 

CO 

C 

o 

co 

CO 

•iH 

CO 

G 

rH 

bD 

O 

O 

XI 

w 

c 

0 

CO 

co 

w 

-P 

•1— 1 

0 

£ 

s 

fa 

a 

T? 

> 

<l> 

o 

H 

CO 

3 

0 

CO 

in 

O 

G 

15 

o 

w 

p 

o 

a 

Eh 

• 

(■  LA 

VO  1 

Propagation  characteristics  @ 75Hz  for  the  segmented  path 
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This  is  a 5.9  dB  signal  strength  reduction  over  ambient  conditions. 

6.2  Calculations  for  Hypothesized  PCA  Profiles 

Using  the  waveguide  computer  program  an  investigation  was  made  of  the 

effects  on  the  ELF  signal  strengths  of  varying  electron  density  profiles 

over  a range  of  values  expected  for  PCA  events.  These  studies  were  performed 

with  variations  of  the  density  profile  based  on  measurements  during  the 

1254  UT  5 August  1972  event  (Reagan  et  al.,  1973).  Since  the  computer  runs 

indicated  that  for  the  densities  of  concern  the  calculated  signal  strengths 

were  dependent  on  the  positive  ion  profile,  for  each  of  the  electron  density 

profiles  the  positive  ion  profile  was  adjusted  so  that  both  are  consistent 

with  the  same  ion-pair  production  profile.  At  all  altitudes  the  ion-ion  re- 

-7  3 

combination  rate  coefficient  was  taken  to  be  10  cm  -sec.  To  check  the 
sensitivity  of  the  calculated  signal  strengths  to  the  positive  ion  densities 
in  the  altitude  regions  of  enhancement  during  PCA  events,  a series  of  runs 

■f 

was  made  for  vastly  differing  N density  profiles  with  a given  electron  den- 
sity profile.  These  are  shown  in  Figure  6-2.  Each  of  the  N+  curves  is 
labelled  with  the  calculated  signal  strength.  Although  an  exhaustive  study 
of  the  sensitivity  of  the  ELF  transmission  to  the  ion  densities  in  each  of 
several  altitude  regions  was  not  performed,  these  runs  suggest  that  the  signal 
strengths  are  most  sensitive  to  the  N+  densities  at  altitudes  of  ~ 45  km  or 
lower.  To  check  the  sensitivity  of  the  calculated  signal  strengths  to  the  ion 

ion  recombination  rate  coefficient,  one  of  the  runs  was  made  with  a coef- 
-8  ^ + 

ficient  of  6 x 10  cm  -sec.  For  that  run  the  N densities  at  altitudes  of 
~ 50  km  were  increased  by  29$,  but  no  changes,  within  0.1  dB,  occurred  in 
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5 AUGUST  1972 


Various  positive  ion  density  profiles  used  in  the  waveguide  mode  calculation 
along  the  electron  density  profile  shown. 


the  signal  levels. 

The  various  electron  and  ion  density  profiles  used  in  the  field  strength 
calculations  are  shown  in  Figure  6-3.  The  electron  density  profiles  differ 
mainly  in  regard  to  the  lower  altitude  to  which  significant  densities  extend. 
The  signal  strengths  calculated  for  each  of  the  profiles  are  listed  in 
Table  6-3.  The  full  set  of  calculations  was  performed  for  a)  the  WTF  trans- 
mitter inside  of  the  FCA  ionization  region  and  b)  the  transmitter  outside 
of  the  region.  The  signal  strengths  were  1.1  to  2.1  dB  lower  when  the  trans- 
mitter was  not  in  the  region  of  enhanced  ionization.  In  order  to  gain  some 
insight  as  to  the  importance  of  certain  features  of  the  density  profiles  the 
calculated  signal  strengths  are  plotted  in  Figure  6-4  as  a function  of  the 

electron  density  at  60  km  and  also  as  a function  of  the  height  of  the 
3 -3 

10  cm  density  value.  The  calculations  clearly  indicate  a decrease  in 
signal  strength  with  an  increasing  electron  density  or  with  a lowering  in 
altitude  of  a given  density.  But  all  of  the  signal  strengths  are  higher  than 
ambient . 

It  appears  that  a major  FCA  event  with  large  fluxes  of  high  energy  protons 
penetrating  to  low  altitudes,  as  on  4 August  1972,  can  cause  significant  at- 
tenuations in  the  ELF  signal  strengths.  However,  during  weaker  PCA  events 
with  lower  fluxes  of  high  energy  protons  the  ELF  signal  strengths  may  be 
either  attenuated  or  enhanced,  depending  upon  the  details  of  the  resulting 
electron  and  positive  ion  density  profiles.  The  variability  of  the  ELF 
transmission  in  response  to  the  shape  and  magnitude  of  the  density  profiles 
is  illustrated  in  Table  (s)  6-4,  6-5,  and  6-6.  Here  the  attenuation  rates 
and  Wait  excitation  factors,  over  each  of  several  segments  on  the  path 


Calculated  Field  Strengths  for  the  Various 
Density  Profiles  based  on  the  PCA  Event 
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Figure  6-4.  For  the  various  PCA  density  profiles  shown  in  Figure  6-2,  the  calculated 
ELF  signal  strengths  are  plotted  as  a function  of  the  altitude  for  an 
electron  density  of  lC>3cm-3  (left  hand  section)  and  as  a function  of  the 
electron  density  at  60  km  (right  hand  section). 
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from  WTF  to  Norway,  vary  in  either  direction  from  ambient  conditions  and  can 
result  in  an  overall  signal  attenuation  or  amplification  depending  upon  the 
details  of  the  propagation  path.  For  comparison  these  propagation  charac- 
teristics are  also  listed  for  various  REP  profiles,  discussed  in  Section  5, 
and  for  sporadic  E layers,  discussed  in  Section  7. 


IJ 


I T- 


6.3  Comparison  of  Effects  from  PCA  and  REP  Events 

The  PCA  event  on  4 August  1972  contained  large  fluxes  of  high  energy 
penetrating  protons  which  produced  intense  ionization  at  altitudes  below 
40  km.  Softer  PCA  events  with  insignificant  fluxes  of  high  energy  protons, 
such  as  observed  on  5 August  1972,  may  produce  enhanced  ionization  only  down 
to  altitudes  of  50-60  km.  In  comparison,  electron  precipitation  events  can 
exhibit  much  higher  ion  pair  production  rates  at  altitudes  above  ~ 60  km, 
but  even  the  most  energetic  REP  events  observed  to  date  do  not  produce  large 
ionization  at  altitudes  much  below  50  km.  The  bremsstrahlung  produced  by  the 
electrons  interacting  with  the  atmosphere  can  penetrate  much  deeper  into  the 
atmosphere  and  can  produce  significant  ionization  down  to  altitudes  below  30 
km.  Comparisons  between  PCA  and  REP  events  in  regard  to  energy  deposition 
are  illustrated  in  Figure  6-5.  The  resulting  electron  densities  are  shown 
in  Figure  6-6. 

We  have  just  seen  that  the  energy  deposition  profiles  for  PCA  and  REP 
events  may  be  similar  but  the  PCA  enhancements  typically  extend  to  lower 
altitudes.  Pronounced  differences  between  these  classes  of  events  are  generally 
associated  with  the  geographic  distribution  and  hence  in  the  contributions 
along  the  various  segments  of  the  transmission  paths.  Polar  cap  absorption 
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production  rates  versus  altitude  for  selected  FCA  and  REP  events 


Electron  density  profiles  resulting  from  the  PCA  and  REP  ion-pair  production  rate 
profiles  shown  in  Figure  6-5. 
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events  typically  cover  the  entire  polar  cap  which  may  or  may  not  include  the 
WTF  transmitter,  whereas  electron  precipitation  events  are  more  confined  to 
the  outer  radiation  belt  region  surrounding  the  polar  cap.  The  ELF  waveguide 
computer  runs  have  shown  that  a more  localized  ionization,  such  as  occurs  in 
REP  events  may  have  less  effect  on  the  transmitted  signal  strengths  than 
does  a widespread  increase  in  the  electron  density.  This  conclusion  was  con- 
firmed by  computing  the  signal  strengths  expected  at  Tromso  with  the  electron 
precipitation  event  at  0401  HT  UT  on  11  March  1975  for  the  increased  ionization 
being  localized  to  a)  the  WTF  transmitter  and  b)  the  entire  path  from  WTF  to 
Tromso.  In  the  former  case  the  signal  strength  was  -157.6  dB  i.e.  the  signal 
strength  was  enhanced  over  ambient  by  1.7  dB.  In  the  latter  case  the  signal 
strength  was  -155-7  dB  compared  to  an  ambient  level  of  -159-3  dB.  A similar 
trend  for  less  signal  enhancement  when  the  REP  effects  were  localized  to 
the  transmitter  was  calculated  for  transmission  to  Greenland. 
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WAVEGUIDE  MODE  CALCULATIONS  FOR  SPORADIC 
E-LAYER  PROFILES 


In  the  previous  study  (imhof  et  al.,  1976)  calculations  with  the  NOSC 

waveguide  computer  program  were  made  to  test  the  sensitivity  of  the  wave 

propagation  to  sporadic  E layers,  E . For  this  purpose  an  E layer  at 

s s 

heights  of  120-125  tan  was  superimposed  on  the  ambient,  nighttime  conditions. 

The  maximum  electron  density  in  the  layer  was  chosen  rather  arbitrarily  as 

was  the  detailed  variation  of  electron  density  with  height.  The  computer 

attenuation  rates  displayed  a strong  frequency  dependence,  reaching  a very 

high  value  of  10.347  dB/Mm  at  80  Hz.  This  very  high  figure  and  the  rather 

narrow  resonance  absorption  peak  between  60  and  85  Hz  gave  indications  that 

the  results  should  be  treated  with  caution.  Accordingly,  it  was  concluded 

that  more  computer  calculations  should  be  made  to  test  the  sensitivity  of 

the  results  to  peculiarities  in  the  electron  density  profile,  the  thickness 

of  the  E layer  as  well  as  its  height  and  maximum  electron  density, 
s 

Based  on  the  aforementioned  conclusions  we  have  made  a series  of  computer 

runs  with  various  hypothesized  sporadic  E-layer  profiles.  Specifically,  the 

height  and  thickness  of  the  E layer  were  varied.  The  attenuation  versus 

s 

frequency  profile  was  found  to  depend  very  strongly  on  the  altitude  of  the 

E layer.  For  a sporadic  E layer  at  ~ 105  km,  as  measured  over  Wallops  Island 
s 

by  Smith  and  Voss  (1977),  &n  attenuation  versus  frequency  curve  was  computed 
at  the  following  frequencies;  45,  75,  80,  85,  90,  95,  100,  110,  120,  130, 

140,  150,  170,  190  and  200  Hz.  The  runs  revealed  no  peaks  in  attenuation, 
perhaps  indicating  their  presence  beyond  the  range  of  frequencies  at  which 
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the  ELF  transmitters  are  operating. 

Instead  of  pursuing  further  the  computer  calculations  with  a sporadic 

layer  at  altitudes  of  ~ 105  km  we  made  a run  with  the  same  E profile  displaced 

s 

upward  to  118  km  and  found  that  the  attenuation  increased  from  O.65  dB/Mm  to 

1.6  dB/Mm  at  75  Hz.  The  sensitivity  to  the  earth's  magnetic  field  strength 

was  investigated  by  changing  B from  0.59^  gauss  to  0.4l  gauss,  whereupon  the 

attenuation  increased  from  1.6  dB/Mm  to  4.02  dB/Mm.  Subsequently  a series  of 

runs  was  made  with  an  E layer  at  an  altitude  of  ~ 122  km,  where  the  original 

s 

runs  were  performed,  and  with  a variety  of  layer  thicknesses,  as  illustrated 
in  Figure  7-1  and  Table  7-1.  For  these  runs  the  earth's  conductivity  at  the 
WTF  antenna  was  taken  to  be  3-2  x 10_l+  MHO/meter  which  was  the  value  used  in 
the  original  runs.  The  attenuation  rates  are  plotted  in  Figure  7-2  as  a 
function  of  frequency  for  the  various  sporadic  E-layer  profiles  shown  in 
Figure  7-4*  Clearly  large  attenuations  are  predicted,  even  for  sporadic  E- 
layers  as  thin  as  ~ 2 km.  The  computer  runs  have  shown  that  the  calculated 
signal  strengths  are  very  sensitive-to  certain  parameters  such  as  altitude 
of  the  Eg  layer  and  the  magnetic  field  strengths.  Since  this  aspect  of  the 
effort  is  not  directly  in  line  with  the  primary  objectives  of  investigating 
the  correlation  between  ELF  propagation  anomalies  and  energetic  particle 
precipitation  as  measured  from  satellites  the  effects  of  Es  layers  were  not 
pursued  further. 
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Section  8 
SUMMARY 


Past  efforts  to  account  for  the  frequent  absorption  anomalies  en- 
countered in  the  nighttime  propagation  of  ELF  radiowaves  have  not  been  con- 
clusive. It  is  realized  that  at  least  some  of  the  variations  encountered  in 
the  absence  of  sunlight  are  due  to  anomalous  iono^>heric  ionization  caased 
by  the  precipitation  of  energetic  electrons  and  protons  into  the  earth's  at- 
mosphere. Recently  under  the  present  contract,  a successful  correlation  was 
made  for  the  first  time  between  direct  particle  precipitation  and  anomalous 
nighttime  ELF  signal  strengths.  However,  a detailed  quantitative  correlation 
between  the  particle  inputs  and  the  ELF  signal  levels  remained  to  be  demon- 
strated. This  report  describes  the  results  of  an  effort  to  perform  quanti- 
tative correlations  between  ELF  anomalies  observed  at  various  stations  around 
the  world  and  precipitating  energetic  particle  inputs  as  measured  from  a 
satellite. 

Calculations  have  been  made  of  the  predicted  signal  strengths  for  trans- 
missions from  the  Wisconsin  Test  Facility  (WTF)  to  each  of  the  receiving 
stations  at  Maryland,  Connecticut,  Greenland,  Norway  and  Italy.  The  greatest 
emphasis  in  these  calculations  has  been  on  the  path  from  WTF  to  Tromso,  Norway. 
This  path  encompasses  high  latitude  regions  of  operational  importance  where 
significant  transmission  anomalies  may  often  result  from  energetic  particle 
precipitation  into  the  atmosphere  during  geomagnetic  storms  and  polar  cap 
absorption  events.  Computations  of  the  ELF  transmission  from  WTF  to  Tromso 
were  performed  using  an  8-segmented  propagation  path.  For  nighttime  ambient 
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conditions  the  calculated  values  at  42  and  75  Hz  were  found  to  be  in  good 
agreement  with  the  values  measured  in  1975  and  1976.  If  the  predicted  values 
are  corrected  for  the  dipole  moment  of  the  N-S  WTF  antenna  using  values 
supplied  by  Dr.  P.  R.  Bannister  then  the  predicted  values  are  slightly  above 
the  measured  averages.  For  the  other  stations  varying  degrees  of  agreement 
were  achieved,  but  since  an  evaluation  of  the  NOSC  computer  program  for 
ambient  conditions  is  not  a primary  objective  of  the  present  effort  the 
matter  was  not  pursued  further.  For  the  present  purposes  it  is  sufficient 
to  have  demonstrated  a general  consistency  between  the  measurements  and  cal- 
culations as  a baseline  for  evaluating  the  effects  of  energetic  particle  pre- 
cipitation. 

For  strong  FCA  events  with  large  fluxes  of  high  energy  penetrating 
protons  pronounced  ELF  signal  strength  attenuations  of  the  order  of  2 to  6 dB 
are  predicted  for  the  path  WTF  to  Trcmso.  However,  during  weaker  FCA  events 
with  lower  fluxes  of  high  energy  protons  the  ELF  signal  strengths  may  be 
either  attenuated  or  enhanced,  depending  upon  the  ion  pair  production  profiles 
and  the  resulting  electron  and  positive  ion  density  profiles.  The  signal 
strengths  were  found  to  be  1.1  to  2.1  dB  lower  when  the  transmitter  was  not 
in  the  region  of  enhanced  ionization. 

The  effects  of  energetic  electron  precipitation  on  nighttime  signal 
transmission  have  been  studied  with  several  different  approaches,  making 
full  use  of  the  satellite  data  and  the  waveguide  computer  program.  An 
initial  correlative  survey  was  performed  with  the  geomagnetic  K index  as  a 
general  indicator  of  particle  precipitation.  For  the  March -April  1976 
period  the  variations  in  the  signal  strengths  at  Norway  and  Maryland  displayed 
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many  similarities  to  variations  in  the  K indices,  with  greater  signal  strengths 
generally  occurring  at  times  of  geomagnetic  disturbance.  However,  this  con- 
sistency may  have  been  misleading  since  the  agreements  for  a time  period  in 
1975  are  much  less  satisfactory  and  at  times  lower  signal  strengths  were  ob- 
served when  the  geomagnetic  activity  was  high.  For  example,  the  ELF  signal 
strengths  at  Maryland  were  below  the  threshold  of  detection  at  times  in  March 
1975  when  intense  relativistic  electron  precipitation  occurred. 

During  the  coordinated  exercises  in  March-April  1976  the  nightly  aver- 
aged signal  strengths  were  found  to  generally  increase  during  periods  when 
the  fluxes  of  precipitating  electrons  > 130  keV  were  higher  on  the  average. 

On  a brood  time  scale  the  various  ELF  receiving  stations  indicated  a fairly 


consistent  pattern  which  may  be  associated  with  the  energetic  electron  pre- 
cipitation. On  a finer  time  scale,  however,  large  inconsistencies  exist  be- 
tween the  various  stations,  suggesting  the  possible  importance  of  other 
effects. 


An  evaluation  of  the  effects  of  REP  events  on  the  ELF  transmission  sig- 
nal strengths  was  made  on  the  basis  of  the  electron  density  profiles  derived 
from  the  measured  spectra  of  precipitating  electrons.  Due  to  the  limited 
spatial  coverage  in  the  satellite  measurements  and  to  the  limitations  of 
the  waveguide  computer  program  in  regard  to  its  inability  to  treat  varia- 
tions of  the  guide  parameters  perpendicular  to  the  propagation  path,  certain 
assumptions  had  to  be  made  regarding  the  spatial  extent  of  the  electron  pre- 
cipitation measured  on  a given  satellite  crossing.  Positive  ion  density  pro- 
files were  obtained  from  the  ion  pair  production  rates  with  an  assumed  ion-ion 

-7  3 -1 

recombination  rate  coefficient  of  10  cm  sec  . For  the  densities  and  altitudes  of 
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concern  the  calculated  signal  strengths  were  found  to  be  significantly  de- 
pendent upon  the  value  of  the  ion-ion  recombination  rate  coefficient  over 
the  range  from  6 x 1D~  to  10~ 'em  -sec”  . The  waveguide  computer  calculations 
were  performed  for  paths  from  WTF  to  each  of  the  receiving  stations  at  Norway, 
Greenland,  Italy,  Connecticut  and  Maryland.  The  results  indicate  that  hard 
electron  precipitation  events  may  cause  a signal  enhancement,  due  mainly  to 
changes  in  the  excitation  factor.  The  consistent  occurrence  of  anomalies  in 
the  ELF  signal  strengths  during  REP  events,  seems  to  be  supported  by  the  co- 
ordinated satellite  ELF  station  measurements.  During  the  intense  events  of 
October-November  1972  lower  ELF  signal  strengths  were  recorded.  Low  signal 
strengths  were  also  associated  with  the  REP  events  on  March  11-12,  1975-  How- 
ever, during  the  weaker  electron  precipitation  events  in  March-April  1976 
there  was  a trend  toward  higher  signal  levels. 

In  summary,  the  following  major  conclusions  have  followed  from  this 
study  effort: 

• From  both  the  satellite  and  ELF  station  measurements  it  has  been 


found  that  direct  particle  precipitation  into  the  atmosphere  can 
cause  ELF  transmission  anomalies.  The  signal  strengths  may  be 


either  attenuated  or  enhanced  depending  upon  the  geometry  and 
details  of  the  ion  pair  production  profiles. 


i 


• Sensitivity  studies  have  been  made  to  assess  the  dependence  of  the 
ELF  signal  strengths  on  such  parameters  as  the  electron  and  ion 
density  profiles  and  their  distribution  along  the  propagation  path. 
The  following  effects  were  discovered: 
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New  techniques  for  mapping  electron  precipitation  profiles 
simultaneously  over  a broad  spatial  region  should  be 
used  in  a coordinated  measurement  program. 

Existing  ELF  waveguide -mode  computer  programs  should  be  modified 
to  include  treatment  of  variations  in  the  electron  and  ion 
density  profiles  along  a direction  perpendicular  to  the  propa- 
gation path. 
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